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ABSTRACT 

A coupling between a scalar field (representing the dark energy) and dark matter could 
produce rich phenomena in cosmology. It affects cosmic structure formation mainly 
through the fifth force, a velocity-dependent force that acts parallel to particle's di- 
rection of motion and proportional to its speed, an effective rescaling of the particle 
masses, and a modified background expansion rate. In many cases these effects en- 
tangle and it is difficult to see which is the dominant one. Here we perform iV-body 
simulations to study their qualitative behaviour and relative importance in affecting 
the key structure formation observables, for a model with exponential scalar field cou- 
pling. We find that the fifth force, a prominent example of the scalar-coupling effects, 
is far less important than the rescaling of particle mass or the modified expansion rate. 
In particular, the rescaling of particle masses is shown to be the key factor leading to 
less concentration of particles in halos than in ACDM, a pattern which is also observed 
in previous independent coupled scalar field simulations. 
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1 INTRODUCTION 

The nature of the dark energy (Copland et al 2006) driving 
an apparent acceleration of the universe has been a cosmo- 
logical puzzle for more than a decade. Models incorporat- 
ing scalar fields are the most popular proposal to explain 
it, those, not only because of their mathematical simplicity 
and phenomenological richness, but also because the scalar 
field is a natural ingredient of many high-energy physics the- 
ories. A scalar field contributes a single dynamical degree 
of freedom which can interact indirectly with other matter 
species through gravity or couple directly to matter, pro- 
ducing a fifth force on the matter which creates violations 
of the Weak Equivalence Principle (WEP). This second pos- 
sibility of direct coupling to matter was introduced with 
the hope that such a coupling could potentially alleviate 
the coincidence problem of dark energy (Amendola 2000) 
and has since then attracted a great deal of attention (see, 
for examples, Bean & Magueijo (2001); Amendola (2004); 
Koivisto (2005); Lee, Liu & Ng (2006); Bochmcr et al. 
(2008); Bean et al. (2008); Bean, Flanagan & Trodden 
(2008); Boehmer et al. (2009) and references therein). It 
was also investigated in the context of theoretical studies 
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of the cosmological variation of the fine structure constant 

(Sandvik, Barrow & Magueijo 2002). 

If there is a direct coupling between the scalar field 
and baryons, then the baryonic particles will experience a 
fifth force, which is severely constrained by observations, 
unless there is some special mechanism to suppress the 
fifth-force effects. This happens in chameleon models, where 
the scalar field (the 'chameleon') gains mass in high-density 
regions (where observations and experiments are performed) 
whereas the fifth force effects are confined to undetectably 
small distances (Khoury & Weltman 2004; Mota & Shaw 
2007). A common approach which avoids such complications 
is to assume that the scalar field couples only to the dark 
matter, an idea seen frequently in models with a coupled 
dark sector [e.g. Caldera-Cabral, Maartens & Schaefer 
(2009); Valiviita, Maartens & Majerotto (2010); 

Simpson, Jackson & Peacock (2010)). In this work our 
scalar field will not be chameleon-like as this case has 
already been investigated elsewhere (Li & Zhao 2009; 
Zhao et al. 2010; Li & Zhao 2010). 

A scalar field coupled to (dark) matter could affect cos- 
mic structure formation in various ways. Firstly, the back- 
ground expansion rate gets modified, which will lead to 
faster or slower clustering of matter particles; secondly, the 
coupling effectively rescales the mass of the particles for the 
coupled matter species, changing the source term of the Pois- 
son equation, which receives a further contribution from den- 
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sity perturbations of the scalar field; thirdly, the coupling to 
the scalar field produces a fifth force between matter parti- 
cles, helping matter to cluster more strongly; finally, there is 
an extra velocity-dependent force on the coupled particles, 
which can be viewed either as a frictional force or as part of 
the fifth force under a frame transformation, and this force, 
being attractive, also promotes matter clustering. 

Not all of these effects are always manifest. In fact, in 
some coupled-scalar-field models one or more then one of 
them could be negligible. An example is the model consid- 
ered in Li & Zhao (2009, 2010), where only the third ef- 
fect, i.e., the fifth force, is non-negligible. This situation 
is also assumed in some other studies of the effect of a 
fifth force on structure formation, e.g., the ReBEL model 
Nusser, Gubser & Peebles (2005), where a Yukawa- type ex- 
tra force is added while the background cosmology is taken 
to be the same as ACDM. 

In other, more general, models, however, a coupling be- 
tween scalar field (dark energy) and (dark) matter often not 
only produces the fifth force, but equally likely creates other 
effects listed above. Examples are the models investigated 
in Maccio et al. (2004); Baldi et al. (2010); Li & Barrow 
(2010), where the modified expansion rate, varying parti- 
cle mass and frictional force are all non- negligible. The sit- 
uation then becomes complicated here, because these could 
have both positive and negative effects on the structure for- 
mation, which are difficult to disentangle. 

In order to clarify the importance of all these effects, we 
have to make detailed analysis by suppressing one or more 
of these effects and then comparing the results - which is 
our primary aim of this work, at least for the model consid- 
ered Li & Barrow (2010). In particular, we would like to see 
how those above effects affect the nonlinear matter power 
spectrum, mass function and profiles of dark matter halos. 
The latter is quite interesting, as it has been shown in ref. 
Baldi et al. (2010); Li & Barrow (2010) that the coupling, 
which is supposed to boost clustering of matter, does in- 
deed suppress the density in the inner region of the halos, 
possibly due to the frictional force and/or the varying mass 
effects. A clarification of the importance of these effects is 
also relevant to the general model tests. As we mentioned 
above, the majority of the investigations of the fifth force to 
date focus only on the fifth force itself but neglect the other 
effects. Yet the same physics responsible for the fifth force 
will often create other associated effects and the latter must 
be taken into account for the sake of consistency. If it turns 
out that the fifth force effects do not dominate the others, 
then a model- independent test of the fifth force might be dif- 
ficult to obtain in cosmology, since different models predict 
very different background expansion and mass variation. 

Since we want to investigate the nonlinear regime, we 
will use the iV-body simulation technique introduced in 
Li & Zhao (2010); Li & Barrow (2010) and also applied in 
Li (2010); Li, Mota & Barrow (2010a, b). Note that other 
approaches to A-body simulations for (coupled) scalar 
field and related models, without solving the scalar field 
equation of motion explicitly, have been also used in var- 
ious previous work, e.g., Maccio et al. (2004); Baldi et al. 
(2010); Linder & Jenkins (2003); Mainini et al. (2003); 
Kesden & Kamionkowski (2006); Springel & Farrar (2007); 
Farrar & Rosen (2007); Keselman, Nusser & Peebles 
(2009, 2010); Hellwing & Juszkiewicz (2009); 



Hellwing, Knollmann & Knebe (2010); Baldi (2010) 
Baldi & Pcttorino (2010); Baldi & Vicl (2010) 
Hellwing, Juszkiewicz & van de Weygaert (2010) 
De Boni et al (2010), while the new feature of our ap- 
proach is that we solve the scalar field equation directly 
(in the quasi-static limit, for more details see Li & Zhao 
(2010); Li & Barrow (2010)). 

The outline of the paper is as follows: in Sect. 2 we 
introduce the basic equations needed to understand the un- 
derlying physics of the model. In Sect. 3 we briefly describe 
the simulations we have performed for the study of the dif- 
ferent effects of a coupled scalar field on structure formation, 
and then display and discuss the numerical results; we sum- 
marise and conclude in Sect. 4. 



2 BASIC EQUATIONS 

All the equations relevant for the simulations used here are 
derived and discussed in detail in Li & Barrow (2010), but to 
make the present work self-contained we list the minimum 
set necessary for us to understand the physical evolution. 
Instead of writing down the field equations directly, we start 
from a Lagrangian 



R 



+ V(ip) - C(<p)C OM + Cs (1) 



where R is the Ricci scalar, k = 8nG with G the gravita- 
tional constant, £dm and Cs ar e respectively the Lagrangian 
densities for dark matter and standard model fields, ip is the 
scalar field, and V(<p) its potential; the coupling function 
C(ip) characterises the coupling between ip and dark mat- 
ter. Given V(ip) and C(ip) a model is then fully specified. 

Varying the total action with respect to the metric 
g a b, we obtain the following expression for the total energy- 
momentum tensor in this model: 



T ab = Va^Vi^ — g a b 



\v c V c y - V(<p) 



(2) 



where T^ 1 and Tf b are the energy-momentum tensors for 
(uncoupled) dark matter and standard model fields. The ex- 
istence of the scalar field and its coupling change the form of 
the energy-momentum tensor, and so modify the cosmology 
from background expansion to structure formation. 

The coupling to scalar field produces a direct interaction 
(a.k.a. the fifth force) between dark matter particles, due to 
the exchange of scalar quanta. This is best illustrated by the 
geodesic equation for dark matter particles 



c(*0 



(3) 



where r is the position vector, t the (physical) time, $ the 
Newtonian potential and V is the spatial derivative. C v = 
dC '/dtp. The second term on the right-hand side is the fifth 
force and only exists for coupled matter species (dark matter 
in our model). The fifth force also changes the clustering 
properties of the dark matter. Note that on very large scales 
ip is homogeneous and the fifth force vanishes. 

In order to solve the two equations above numerically 
we need to solve both the time evolution and the spatial 
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distribution of p, and this could be done using the scalar 
field equation of motion 



„ dVjip) dC{ip) 
V V a ip H ; h Pdm — ; = U 

dip dtp 



or equivalently 



dtp 



where we have defined 



Veff(<fi) = V(<p) +PdmC(v?). 



(4) 



(5) 



(6) 



The background evolution of ip can be solved easily once 
we know the current pdm, because pdm oc a~ 3 . We can 
then divide ip into two parts, ip = ip + Sip, where <p is the 
background value, and 5<p the (not necessarily small and 
linear) perturbation, and subtract the background scalar- 
field equation of motion from the full equation to obtain the 
equation of motion for Sip. In the quasi-static limit where we 
can neglect time derivatives of Sip compared with its spatial 
derivatives (which turns out to be a good approximation for 
our simulations, because the simulation box is much smaller 
than the observable Universe), we get 

dC(p) dV(ip) dV(p) 

jz — Pdm H j — — (7) 



V ip = — — — pdm 



dip dip dip dip 

where Pdm is the background dark-matter density. 

Once Pdm is known on a grid, we can then solve 
Sip on that grid using a nonlinear Gauss-Seidel relax- 
ation method (in our simulations we have modified MLAPM 
Knebe, Green & Binney (2001), a publicly available TV-body 
code using a self-adaptive refined grid so that high resolu- 
tions can be achieved in high-density regions). Since <p is also 
known, we can then obtain the full solution of ip = ip + Sip. 
This completes the computation of the source term for the 
Poisson equation: 

V 2 $ = | [<7(<p)p D M - C{ip)p DM + 5p B - 25V(p)] , (8) 

where SpB = Pb— Pb and SV(ip) = V(ip) — V(<p) are respec- 
tively the density perturbations of baryons and scalar field 
(note that we have neglected perturbations in the kinetic 
energy of the scalar field because they are always very small 
for our model). 

We can then solve Eq. (8) using a linear Gauss-Seidel 
relaxation method on the same grid to obtain $. With both 
$ and ip in hand, Eq. (3) can then be used to compute the 
forces on the dark matter particles, and once we have the 
forces, we can perform all the standard TV-body operations 
such as momentum-kick, position-drift, time-stepping and 
so on. 

Eqs. (2 - 8) are all what we need to complete an TV-body 
simulation for coupled scalar field cosmology Li & Barrow 
(2010), and from them we can see where the effects of the 
scalar-coupling enter: 

(i) The modified background expansion rate mainly af- 
fect the particle movements and time-stepping, i.e., Eq. (3), 
because in the simulations we use the cosmic scale factor a, 
instead of t, as the time variable and d/dt = ad /da. 

(ii) The varying mass effect is seen directly from Eq. (8), 
which shows that the contribution of pdm to the source term 
of the Poisson equation is normalised by C(ip) which is dif- 
ferent from 1 in general. In our model it is not true that the 



mass of dark matter particles is really varying, but the net 
effect is just equivalent to such a variation. 

(iii) The fifth force appears explicitly on the right-hand 
side of Eq. (3), but is only for coupled matter species (dark 
matter). 

(iv) The velocity-dependent (or frictional) force hides in 
the fact that Eqs. (3, 7) are given in different gauges: Eq. (3) 
is the force for a dark-matter particle and is given in that 
particle's rest frame, while Eq. (7) is written in the funda- 
mental observer's frame. As a result, in order to use the Sip 
solved from Eq. (7) in Eq. (3), we need to perform a frame 
transform VSp — > VSip + aip-k, where x is the comoving ve- 
locity of the particle relative to the fundamental observer. 
This force is therefore expressed as —^r-aipji, and obviously 
the faster a particle travels the stronger such force it feels. 



3 SIMULATIONS AND RESULTS 

3.1 Model and Simulation Details 

As is mentioned above, a coupled scalar field model is fully 
specified given the exact forms of the bare potential V(tp) 
and coupling function C(<p). Here we shall choose the same 
model as Li & Barrow (2010), with an inverse power-law 
potential 



V(ip) 

and the exponential coupling 



A 4 



(9) 



C(p) = exp (jy^p) , (10) 

where A is a constant with mass dimension and A 4 is of 
order the dark-energy density today; a, 7 are dimensionless 
parameters. We choose a — 0.1 so that the potential is flat 
enough to enable a slow-roll of ip (which accounts for the 
dark energy); for 7, we choose I7I ~ O(0. 1) and 7 < so 
that both V(ip) and V e //(<p) are of runaway type 1 . 

Li & Barrow (2010) have given a very detailed descrip- 
tion of the technicalities of the TV-body simulations for cou- 
pled scalar-field theories, and so we shall not repeat this 
methodology here. Roughly speaking, the most important 
distinction between our simulation and others is that we 
have solved the scalar field equation of motion explicitly on 
a grid (in the quasi-static limit). Because of this, we have 
been able to solve the fifth force (and the frictional force) 
numerically without recourse to analytical approximations. 
Furthermore, we have incorporated both time and the space 
variations of the particle mass (or more rigorously of C(ip)) 
because we have spatial information about the scalar field 
distribution. 

Since our aim is to test the significance of each of the 
four above-named effects, we choose to suppress one of them 
at one time. Together with the full model, where all effects 
are included, and the ACDM model for comparison, we then 
have six models to simulate. Furthermore, we consider two 



1 Since V e ff(tp) is of runway type there is nothing to stop 
the scalar field rolling down V e f f , so typically we shall have 
i/Rip ~ O(l) today, which makes C(ip) deviate significantly from 
1. Increasing I7I will makes this problem more severe, and this is 
why we set \y\ ~ O(0. 1) rather than 0(1). 
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simulation no. 


a 


7 


simulation description 


L 


0.0 


0.00 


pure ACDM 


SI 


0.1 


-0.10 


full coupled scalar field 


Sla 


0.1 


-0.10 


scalar field with the frictional force — ^-a(pii suppressed 
scalar field with the fifth force — ^-\/5ip suppressed 


Sib 


0.1 


-0.10 


Sic 


0.1 


-0.10 


scalar field with the (time and spatial) variation of mass C(ip) removed 


Sid 


0.1 


-0.10 


scalar field with a ACDM background expansion 


S2 


0.1 


-0.20 


full coupled scalar field 


S2a 


0.1 


-0.20 


Q 

scalar field with the frictional force — ^-ai^x suppressed 


S2b 


0.1 


-0.20 


C 

scalar field with the fifth force — j#-V<5i/j suppressed 


S2c 


0.1 


-0.20 


scalar field with the (time and spatial) variation of mass C(ip) removed 


S2d 


0.1 


-0.20 


scalar field with a ACDM background expansion 



different choices of the coupling strength 7. So in total we 
have 11 models, details of which are summarised in the fol- 
lowing table: 

The physical parameters we adopt in all simulations 
are as follows: the present-day dark-energy fractional energy 
density Q DE = 0.743 and D. m = S1 C dm + Ob = 0.257, H = 
71.9 km/s/Mpc, n 3 = 0.963, a 8 = 0.761. The size of simula- 
tion box is 64/i _1 Mpc with h = H o /(W0 km/s/Mpc). In all 
these simulations, the mass resolution is 1.114x 10 9 /i _1 Mg, 
the particle (both dark matter and baryons) number is 256 3 , 
the domain grid is a 128 x 128 x 128 cubic and the finest 
refined grids have 16384 cells on each side, corresponding to 
a force resolution of order 12/i _1 kpc. 

3.2 Numerical Results 

For the numerical results, we display the matter power spec- 
trum, mass function and halo density profiles for the 11 runs 
described above, and discuss how they are affected by the 
individual effects from the scalar coupling. 

Before going to the details, it is helpful to have a quick 
browse about the scalar-coupling effects: (I) The velocity- 
dependent force — ^-aifiX is parallel to the direction of mo- 
tion (x) ; because <p > and 7 < 0, so it accelerates the par- 
ticles. (II) The fifth force — ^rVStp in this model is found 
(Li & Barrow 2010) to be parallel to gravity and the ratio 
between the magnitudes of the two is 27 2 to a high preci- 
sion; as such the fifth force both accelerates the particles and 
increases their mutual attraction. (Ill) C(tp) < 1 because 
7 < and <p > 0, so that the contribution of dark matter 
density to the source of Poisson equation gets weakened, ef- 
fectively reducing the gravity force and causing less mutual 
attraction between particles and less clustering. (IV) For the 
chosen model and physical parameters, the background ex- 
pansion rate decreases as I7I increases (Li & Barrow 2010), 
making matter particles less diluted and cluster more. These 
facts are important to bear in mind for discussions below. 

3.2.1 Matter Power Spectrum, 

In Fig. 1 we have plotted the fractional change of the nonlin- 
ear matter power spectrum with respect to ACDM predic- 
tion for the simulations SI and Sla-d. Roughly, the deviation 
from the black solid curve (SI) indicates the importance of 
a given coupled-scalar effect: the larger the deviation is, the 
more that specific effect contributes to the full coupled scalar 



field result. Of least importance is the velocity-dependent 
force —^-a(p'k (green dotted curves). As this force acceler- 
ates particles, it makes particles collapse faster to the regions 
of high density, and thus slightly enhances the clustering. 
Consequently, suppressing it will decrease P(k). Of the sec- 
ond least importance is the fifth force term — ^-VScp (blue 
dashed curves). It not only accelerates particles (towards the 
high-density regions) but also increases the central force. As 
such its effect is also to enhance the clustering of matter and 
neglecting it leads to smaller P(k) on all scales. 

Next is the varying mass effect (purple dash-dot curves), 
the presence of which reduces the source of the Poisson equa- 
tion and thus weakens gravity. Obviously its effect is to pro- 
duce weaker clustering of matter particles and dropping it 
will increase P(k) significantly. 

he most important coupled-scalar effect comes from the 
modified background expansion rate (pink dash-dot-dot-dot 
curves). As mentioned above, if the universe expands more 
slowly (as in the case of our coupled scalar field model 
Li & Barrow (2010)), then particles are less diluted and have 
more time to cluster, resulting in a larger P(k). Changing 
the background expansion to ACDM (i.e., increasing it) sim- 
ply produces a smaller P(k) than the full simulation. Note 
that the pink curves are consistently below zero, indicating 
that although simulation Sid uses the same background ex- 
pansion rate as simulation L, P(k) is smaller than for the 
latter, a result that is again due to the fact that in Sid 
the varying mass effect is taken into account, weakening the 
matter clustering and decreasing P(k) (as discussed above, 
the fifth force does have the opposite effect but cannot over- 
come this). 

In Fig. 2 we have shown the same plots, but for the 
models S2 and S2a-d. All the above analysis still applies but 
the effects just become stronger. Interestingly, the most im- 
portant coupled-scalar effects (at least for our models, which 
are typical ones) are not those of the fifth force, but are in- 
stead the modified background expansion rate and varying 
mass of matter particles. This is understandable, because 
the magnitude of the fifth force is about 27 2 times that of 
gravity, and for our SI and S2 simulations the values are 
0.02 and 0.08 respectively, while at the same time the devi- 
ations of C((f) from 1 for these two models are ~ 0.1 and 
~ 0.3 Li & Barrow (2010) and so are far stronger. 

This means that one must be cautious about adding 
a Yukawa-type fifth force to the A^-body simulation, while 
keeping all other things the same as in ACDM, because the 
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Figure 1. (Colour Online) Fractional changes of the nonlinear matter power spectrum with respect to the ACDM result at four different 
output times: a = 0.3 (upper left panel), a = 0.5 (upper right panel), a = 0.7 (lower left panel) and a = 1.0 (lower right panel). In each 
panel the results of the simulations SI (full coupled scalar simulation), Sla (frictional force suppressed), Sib (fifth force suppressed), Sic 
(mass variation removed) and Sid (ACDM background) arc represented by the black solid, green dotted, blue dashed, purple dash-dot 
and pink dash-dot-dot-dot curves respectively. 



fifth force often introduces associated effects which are more 
important 



3.2.2 Mass Functions 

Fig. 3 shows the mass functions of the simulations SI and SI 
compared with that of ACDM (L). Again, the deviation from 
the full simulation result (black solid curve) indicates the or- 
der of importance of the individual effects. The two least im- 
portant factors are once more the velocity-dependent force 
— jj-aip'k and the fifth force — jj-X?8ip, both of which, ac- 
cording to our above analysis, enhance matter cluster: sup- 
pressing them causes less clustering of matter and smaller 
mass functions. Their influences however are quite weak, in 
particular that of the velocity-dependent force. 

The second important effect is the variation of mass 
or C(ip). In case C(ifi) < 1 then, as mentioned above, the 
mutual gravitational interaction between particles becomes 
weaker and particles will cluster less strongly. As a result, re- 



moving this effect enhances the matter clustering and leads 
to massive halos being created in larger abundance. 

The most influential effect from a coupled scalar field is 
the modified background expansion. Changing it to a ACDM 
background (which is faster) significantly underestimates 
the mass function, because particles get more diluted and 
have less time to clump. 

The corresponding results for models S2 and S2a-d are 
summarized in Fig. 4, and they show the same qualitative 
trend as Fig. 3 but the effects are just stronger, due to the 
stronger coupling |7| and thus more dramatic evolution of 
the scalar field <p (which means that C(<p) deviate more from 
unity and the expansion rate is decreased more compared to 
simulation L). 

Although the mass function does show the imprint from 
a fifth force, this is by no means unique and could easily 
be dominated over by the associated coupled scalar effects 
like modified background expansion rate or varying particle 
mass. 



6 Baojiu Li, John D. Barrow 





Figure 2. (Colour Online) The same as Fig. 1, but for the models where 7 = —0.20, i.e., S2 and S2a-d. 



3.2.3 Halo Density Profiles 

Internal density profiles of the dark matter halos are another 
area where the scalar field coupling could leave interesting 
imprints. For example, Li & Zhao (2010) showed that for a 
chameleon-like scalar field model the density profiles could 
be either increased or decreased significantly by the scalar 
coupling, depending on the environment of the halos. For 
the models considered here, Baldi et al. (2010); Li & Barrow 
(2010) have given convincing evidence that the internal den- 
sity profile has similar shape to ACDM, but could be some- 
what suppressed in the very inner parts. 

As an explicit example, Fig. 5 shows the density profiles 
of the most massive halo from each box for simulations L, 
SI and Sla-d. The suppression of the full simulation result 
(black solid curve) compared with the ACDM prediction 
(thick solid dashed curve) is evident below R ~ 400/i _1 Kpc. 

If the velocity-dependent force — ^r-atpx is given up 
(green dotted curve), then the suppression is moderated. As 
discussed in Baldi et al. (2010), this is because this force ef- 
fectively accelerates particles, making them travel faster and 
so increases the total kinetic energy of the particles inside 
the halo. As a result, removing this force will leave the parti- 
cles with less kinetic energy, meaning that they more easily 



fall towards the halo centre. Unlike in Baldi et al. (2010), 
here this effect is not the major one (see below). 

The fifth force — -^-V8(p is slightly more complicated. 
On one hand, it speeds the particles up, increasing their ki- 
netic energy; on the other, it enhances the mutual attraction 
between particles, increasing the (magnitude of the) poten- 
tial energy of the halo. These two effects are just opposite, 
as deepening the potential of a halo will pull more particles 
towards the centre and thus increase the concentration. For 
simulation SI, it seems that the former effect takes over, 
and dropping the fifth force (blue dashed curve) simply de- 
creases the kinetic energy of particles and make them more 
concentrated towards the halo centre. 

The effect of the varying-mass factor C(ip) < 1 is un- 
ambiguously to suppress the over-density inside the halos, 
as it decreases the source of the Poisson equation and so 
weakens gravity, making the potential shallower. Dropping 
it simply deepens the potential and attracts more particles 
to the inner region of the halo (purple dash-dotted curve). 

Finally, changing (increasing) the background expan- 
sion rate relative to ACDM again leads to less clustering of 
particles, and thus lower density profiles in the halos (pink 
dash-dot-dot-dot curves) . 
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Figure 3. (Colour Online) The mass functions for the models SI 
(black solid curve), Sla (green dotted), Sib (blue dashed), Sic 
(purple dash-dot) and Sid (pink dash-dot-dot-dot) compared to 
the ACDM result (thick long dashed curve). The horizontal axis 
is the virial mass of halos, in unit of /i _1 AfQ; the vertical axis 
is the halo number density in the simulation box, in units of 
(/i^Mpc) -3 . 
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Figure 5. (Colour Online) The internal density profile of the 
most massive halo in the simulation box for the models SI (black 
solid curve), Sla (green dotted), Sib (blue dashed), Sic (pur- 
ple dash-dot) and Sid (pink dash-dot-dot-dot) compared to the 
ACDM result (thick long dashed curve). The horizontal axis is 
the radius from halo centre, in unit of h~ 1 Kpc, and the vertical 
axis is the overdensity. 
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Figure 4. (Colour Online) The same as Fig. 3, but for models 
S2 and S2a-d. 

It is worth noting that although the velocity-dependent 
force, fifth force, and varying mass all contribute to suppress- 
ing the internal density of the halo, the varying mass con- 
tributes most, while the fifth force contributes least. This dif- 
fers from the finding of Baldi et al. (2010), that the velocity- 
dependence force is the determining factor, possibly due to 
differences in the models and their treatments. The fact that 
the velocity-dependent force dominates over the fifth force 
(which is different from what we have seen for matter power 
spectrum and mass function) is not surprising, for, as men- 
tioned above, the fifth force is set to increase both the kinetic 
and potential energies, two effects somehow cancelling each 
other out. 

Fig. 6 displays the same results for the models S2 and 
S2a-d, and we can see qualitatively similar but stronger 
trends. The notable thing is that in this case suppressing 



10 s 



•■ I I 

- \ ^ 


I I I I I I I I I I I L 




0t=0.1 , 7 = 0.20 










































• 111 


i ivf 



10 100 1000 

R (h" kpc) 

Figure 6. (Colour Online) The same as Fig. 5, but for models 
S2 and S2a-d. 

the fifth force further lowers the inner density of the halo, 
an indication that here the deepening of potential dominates 
over the increase in kinetic energy (both are due to the fifth 
force) . 

Above we have just discussed the results for one specific 
halo (the most massive one), while what we are more inter- 
ested in is the general behaviour. For this we have selected 
10 out of the most massive halos and computed their aver- 
age density profile. The results for simulations SI and Sla-d 
are shown in Fig. 7, in which we have plotted the fractional 
change of the averaged halo internal density with respect to 
the ACDM result. From this figure we see again that overall 
the internal density in the inner region of halos is lower in 
(full) coupled scalar field models than in ACDM. Both the 
velocity-dependent force and varying mass tend to decrease 
the inner density, and both the fifth force and the modified 
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Figure 7. (Colour Online) The fractional change of averaged 
halo density profile for the models SI (black solid curve), Sla 
(green dotted), Sib (blue dashed), Sic (purple dash-dot) and 
Sid (pink dash-dot-dot-dot) compared to the ACDM result (thick 
long dashed curve). The horizontal axis is the distance from halo 
centre, in unit of h~ x Kpc, and the vertical axis is the fractional 
change of density with respect to the ACDM prediction. 




Figure 8. (Colour Online) The same as Fig. 7, but for models 
S2 and S2a-d. 



(slower) background expansion help increase it. The effects 
from the varying mass and modified background expansion 
are dominant while the other two effects are minor, confirm- 
ing our observed pattern from a single halo (Fig. 5). 

We see a similar result for the models S2 and S2a-d, as 
shown in Fig. 8. A notable difference here, however, is that 
the velocity-dependent force becomes more important than 
in SI. Indeed, it is as dominant as the varying mass effect in 
the inner region of the halos, though it is still subdominant 
in the outer region; this at least shows some agreement with 
Baldi et al. (2010). Note again that the fifth force is the 
least important amongst all the four effects, which means 
that using a simulation with ACDM plus a Yukawa-type 
fifth force we would be unable to catch the (probably) most 
significant effects from a coupled scalar field. 



To summarise: in this paper, with the aid of iV-body sim- 
ulations, we have investigated the different impacts of a 
coupled scalar field on cosmic structure formation, and as- 
sessed their qualitative effects and quantitative importance. 
The scalar field coupling influences the structure formation 
mainly through a velocity-dependent force, a fifth force, a 
modification of the particle mass (or the source of the Pois- 
son equation) and a modified background expansion rate. 
We have investigated how dropping each one of these factors 
leaves imprints on the key structure formation observables, 
like matter power spectrum, mass function and the internal 
density profile of dark matter halos. 

For the matter power spectrum and mass function, we 
find that the modified background expansion rate is by far 
the most important effect that the scalar-coupling can have, 
followed in turn by the variation of particle mass, fifth force 
and the velocity-dependent force. The cosmic expansion be- 
comes slower than that in ACDM due to the scalar field, 
which means that structure has more time to form. The 
fifth force increases the mutual attractive force between par- 
ticles, strengthening the collapse of overdense regions; both 
the fifth force and the velocity-dependent force could speed 
up particles, making the collapse faster. As a result, all these 
three effects help boost the growth of structure. On the other 
hand, the source of the Poisson equation is decreased due to 
the coupling function C(ip), resulting in weaker gravity and 
weakened structure formation. 

The internal density profiles for dark matter halos are 
more interesting, and as we have seen the combined effect of 
a scalar-coupling can be to suppress the density of the inner 
regions of the halos (or to distribute particles more towards 
outer regions). As one might expect, the modified cosmic 
expansion rate is again the most important scalar-coupling 
effect here, and with it dropped there will be less structure 
formation and lower halo density profiles. The variation of 
particle masses, or the modification to the source of the Pois- 
son equation, is the second largest single effect, followed by 
the velocity-dependent force and fifth force. Roughly speak- 
ing, the particles tend to move towards (by a process of re- 
laxation and virialisation) the inner regions of halos if their 
kinetic energy is reduced and/or the central potential gets 
deeper, and vice versa (Baldi et al. 2010). In this regard, the 
velocity-dependent force speeds up particles and the varia- 
tion of particle mass weakens the central potential, both in 
favour of lower central densities in halos (cf. Figs. 5-8). The 
fifth force has two opposite effects - to speed up particles 
(and so increase the kinetic energy) and to deepen the total 
potential. They cancel each other out to a certain extent, 
so making the fifth force less influential in determining the 
density profiles (than the velocity-dependent force) 

Our result is marginally consistent with previous analy- 
ses (e.g., Baldi et al. (2010)) of the halo density profiles, but 
shows discrepancy about whether the velocity-dependent 
force is more important than the varying particle mass or 
not. The difference might be due to different model specifi- 
cations and treatments. 

One of our most important results is that in many cases 
the fifth force, which is the most well-known consequence 
of a coupling between matter and a scalar field, is not the 
most important in affecting the structure formation. The 
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key point is that, when we introduce such a coupling, other 
new effects are also brought in, and these can often be much 
more influential. It is in this regard that the advantages 
of full iV-body simulations Maccio et al. (2004); Baldi et al. 
(2010); Li & Zhao (2009, 2010); Li & Barrow (2010), which 
take full account of all associated effects, are increasingly 
significant. 
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